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Abstract 

Based upon elements of the modern Pseudoanalytic Function Theory, 
we analyse a new method for numerically approaching the solution of the 
Dirichlet boundary value problem, corresponding to the two-dimensional 
Electrical Impedance Equation. The analysis is performed by interpolat- 
ing piecewise separable- variables conductivity functions, that are eventu- 
ally used in the numerical calculations in order to obtain finite sets of 
orthonormal functions, whose linear combinations succeed to approach 
the imposed boundary conditions. To warrant the effectiveness of the 
numerical method, we study six different examples of conductivity. The 
boundary condition for every case is selected considering one exact solu- 
tion of the Electrical Impedance Equation. The work intends to discuss 
the contributions of these results into the field of the Electrical Impedance 
Tomography. 

1 Introduction 

The study of the Dirichlet boundary value problem for the two-dimensional 
Electrical Impedance Equation 

div(crgradu) — 0, (1) 

is fundamental for well understanding its inverse problem, usually know as Elec- 
trical Impedance Tomography, and that was first posed in correct mathematical 
form by A. P. Calderon [3] in 1980. It is remarkable that for more than twenty 
years after the problem was stated, the mathematical complexity of ([T]) could 
provoke that many experts considered impossible to obtain its general solution 



in analytic form [TH], even for the simplest cases of a. But the this conception 
turned around in 2005, when V. Kravchenko [TT] noticed that the two dimen- 
sional case of ([T]) was completely equivalent to a special kind of Vekua equation 
|T7j. Only one year latter, K. Astala and L. Paiviirinta [T] rediscovered indepen- 
dently this relation, and with it, they gave a positive answer to the Calderon's 
problem in the plain. 

Many other important results were obtained soon after. As an example, 
Kravchenko et al. published in 2007 what can be considered the first general 
solution of ([T]) in analytic form [12] , proposed for a certain class of conductivities 
(T, and employing Taylor series in formal powers [5]. But we would like to claim 
now special attention on a recent and relevant discovering: The proof of the 
completeness of the set of formal powers in bounded domains, for approaching 
solutions in the plain of the Dirichlet problem for the Electrical Impedance 
Equation |5j . This will be valid when the conductivity function can be expressed 
by means of a separable- variables function. 

On the light of all these notorious advances, we want to emphasize that it is 
not clear yet how to fully apply the elements of the modern Theory of Pseudo- 
analytic Functions [5] into more specific, but still very important. Engineering 
applications, as it is the Electrical Impedance Tomography from the point of 
view of Medical Imaging. 

The main objective of this work is to start a discussion in this direction. 
Even it is not clear if the proof provided in |5j can be extended for the cases 
when CT is a separable- variables function, but piecewise-defined within a bounded 
domain, the numerical calculations show that the techniques employed in purely 
mathematical problems, could well be used for analysing some physical cases, 
providing quite acceptable results. 

Basically, our work suggest that, if the values of the electrical conductivity 
are known at every point within a bounded domain, it will be always possible to 
introduce a piecewise separable-variables function that approaches the conduc- 
tivity, and in consequence, to obtain a numerical orthonormal set of functions 
that seems to be capable of approaching solutions for the Dirichlet boundary 
value problem of ([T]). This would be true for a certain class of bounded do- 
mains defined on the plain, but wide enough to include most relevant cases for 
Experimental Physics. 

We base our assessments onto a set of conductivity examples that tries to 
be as diverse as possible from the mathematical point of view. Possessing for 
each example an analytic solution, we introduce their corresponding piecewise 
separable- variables conductivity functions, and employing them for the numer- 
ical calculations, we approach their corresponding analytic solutions valued at 
the boundary. Then, a variety of qualitative and quantitative references are 
given, in order to estimate the effectiveness of the approach. 

We make special emphasis in those examples whose conductivity functions 
are not separable- variables by definition, showing that, at least for the analysed 
cases, the technique can be trustful. This would be an important step for 
performing experiments with other classes of conductivities, as those emerging 
from geometrical distributions, whose exact mathematical expressions are, in 
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general, unknown. If our appreciations are correct, we would be one step closer 
for fully applying the modern Pseudoanalytic Function Theory into Medical 
Imaging. 



2 Preliminaries 

Following [2], let the pair of complex- valued functions F and G fulfil the condi- 
tion 

Im(FG) > 0, (2) 

where F represents the complex conjugation of -F: F ~ Rei^ — zlmi^, and i is 
the standard imaginary unit = —1. Therefore, any complex- valued function 
W can be written by means of the linear combination of the pair (i^, G): 

W = (I)F + tPG, 

where and tp are purely real-valued functions. Based upon this idea, L. Bers 
introduced the concept of the (F, G)-derivative of W as 

d(F,G)W^{d,<l>)F+{d,^)G. (3) 

But this derivative will exist if and only if the following condition holds: 

(chcb) F + G = 0. (4) 

From hereafter, we will consider dz ~ dx~idy and dz = dx+idy. Nevertheless 
these operators are usually introduced with the factor ^, in this work will be 
somehow more convenient to work without it. 

By introducing the notations 

_FdzG~GdzF _ FckG-GdzF 

V.G)- ^^-GF ' FG-GF ' 

„ FdzG-GdzF Gd^F-FdzG 

FG~GF ' ^^^G)-- ; (5) 

the (F, G)-derivative of W, presented in ([3|, can be written as 

diF,G)W = dzW-A^F^G)W-B^F^G)W, (6) 

whereas the condition Q will turn into 

dzW-a^F,G)W-b^F,G)W = 0. (7) 

A pair of complex functions (F, G) satisfying ([2|, will be called a generating 
pair, and the notations introduced in ^ will be referred as the characteristic 
coefficients of the generating pair (F, G). Moreover, the expression ([t]) is known 
as the Vekua equation |17) . and soon it will be seen that it possesses special 
significance for this work. It is also necessary to mention that every function 
W, solution of ([t]), will be called {F,G) -pseudoanalytic. 

The following statements were originally posed in [2] . We present them here 
slightly modified in behalf of better explaining our results. 
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Theorem 1 The elements F and G, belonging to the generating pair are 
both {F, G) -pseudoanalytic and their {F,G) -derivatives ^ vanish identi- 
cally 

d{F.,G)F = d[F,G)G = 0. 

Theorem 2 JE/l^ Let p be a non-vanishing function within some domain f2, 
and let 

F^p, G=-. (8) 
P 

It is easy to verify that F and G constitute a generating pair whose char- 
acteristic coefficients ^ are 

^(F,G) = «(-F,G) — 0, 
^{F,G) = —, O^F,G) = ~- 

Therefore, the corresponding Vekua equation ^ of an {F, G) -pseudoanalytic 
function W will have the form 

d^W - —W = 0. (10) 

P 

Definition 1 Let (_fo,Go) and {Fi,Gi) be two generating pairs of the form 
and let their characteristic coefficients fulfil the condition 

B{Fo,Go) -&(Fi,Gi)- 

The generating pair (Fi, Gi) will be then called a successor pair of {Fq, Gq), as 
well (Fq, Go) will be named a predecessor pair of {Fi, Gi). 

Definition 2 Let the elements of the set 

{{F^,Gra)}; m = 0,±l,±2,... (11) 

be all generating pairs, and let every {F„i^i,G„i+i) be a successor {F„i,Gm)- 
Hence, the set (11 ) will be called a generating sequence. Particularly, if {F, G) = 
{Fq, Gq), we will say that {F, G) is embedded into ( 11). Moreover, if there exist a 
number k such that (Fm+k, Gm+k) — (-Fm, Gm), we. will say that the generating 
sequence (11) is periodic, with period k. 

L. Bers also introduced the concept of the [F, G)-integral of a complex func- 
tion W . The complete and detailed conditions for its existence can be found in 
P] and in [9], but they are out of the scope of this work. Because of this, the 
elements we present in the upcoming paragraphs are only basic statements for 
supporting our results. 

Definition 3 Let {Fq,Gq) be a generating pair with the form Its adjoint 
pair {Fq,Gq) will be defined as 

f; = -iFo, G; = -iG„. 
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Definition 4 The {Fq,Gq) -integral of a complex function W (if it exists) is 
defined according to the expression 

I Wd(Fo,Go)Z^FoRe I GlWdz + G^Re [ F*Wdz, 

Ja ' J A J A 

where A is a rectifiable curve going from Zq upto zi , in the complex plain. In 
particular, the {Fq^Gq) -integral of d(^pgCo)^ reaches 

d(^Fo,Go)WdF,^GoZ = W- <j>{zo)Fo - ijizo)G„. (12) 



But according to the Theorem [7| the {F^^^G^)- derivatives of Fq and Gq vanish 
identically, hence the equation pM) can be considered the {Fq, Gq) -antiderivative 

ofd(Fo,Go)W. 

2.1 Formal Powers 

Definition 5 The formal (a„, Zq] z) belonging to the generating pair 

{Fm,Gm), with formal exponent n, complex coefficient an, center at zq, and 
depending upon the complex variable z = x -\- iy, is defined by the expression 

Z^f> (a„, zo;z) = \F„i + tiG,n, 

where A and fi are complex constants that fulfil the equality 

XFmizo) + ^iGrnizo) = ao- 

The formal powers with higher formal exponents, are defined according to the 
recursive formulas 

Z^^ (a„, Zo; z)^n I Z^^Zp (a™, zq; z) d(F„,,G,^)Z- (13) 

Notice the integral operators at the right hand side of the last equality, are all 
{F„i, Gm)-antiderivatives. 

Remark 1 The formal powers posses the following properties: 

1. zln'' (an, Zo; z) a„ (z - Zo)" when z zq. 

2. All Zm"^ (an, Zq; z) are (F„i,Gm)-pseudoanalytic. 

3. If a^ — a'^ -\- ia'^, where and a^ are real constants, we will have that: 

Z^m (an, Zo; z) = a^Z,(") (1, zo; z) + a"Z^''> {i, zq; z) . (14) 
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Theorem 3 Every complex-valued function W , solution of the Vekua equation 
Q), can be expanded in therms of the commonly called Taylor series in formal 
powers: 



(15) 



where the absence of the subindex "m" indicates that all formal powers belong 
to the same generating pair. 

Remark 2 Since every W , solution of accepts the representation ( 15), it is 
possible to assert that {15) is an analytic representation of the general solution 
for the Vekua equation^T^^- 



3 The Electrical Impedance Equation 

As it has been previously posed in several works (see e.g. [S], [3], [H] and 
|13j). when the conductivity function a in ([T]) can be expressed by means of a 
separable-variables function 

cr = ai{x)a2{y), (16) 

when introducing the notations 

W = \fadr^u — i^/adyU, 



P=^; (17) 



the two-dimensional Electrical Impedance Equation ([T]) can be rewritten pre- 
cisely as a Vekua equation of the form ( 10 ). Moreover, its corresponding gener- 
ating pair 

Fq Go = -, 

P 

is embedded into a periodic generating sequence, with period 2, such that 
1. 



F G -ly""^ 



when m is an even number, and 

2- 

^ 

Fm — -v/oTa/o^, G„i = 



when m is odd. 
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Therefore, on the hght of the statements posed in Definition |5] the exphcit 
generating sequence will allow us to approach the set of formal powers: 

(4"\l,0;z),z(")(z,0;z)r , (18) 

^ J n— 

within a bounded domain Q,, and by virtue of the Remark [T] we will be able 
to approach any formal power ZQ"^(a„, 0; z); n — 0, 1,..., N] at some certain 
point z € Q. 

Since the present work intends to be part of a novel theory for the Electrical 
Impedance Tomography problem, we will focus our attention into a classic do- 
main f2: The unitary disk with center at zg = 0- We shall point out that all the 
results we will present, can be generalized for a wide class of bounded domains. 
However, it will be more comfortable by now to perform our analysis into this 
classical domain. 

Hereafter, we will analyse the Dirichlet boundary value problem for the two- 
dimensional Electrical Impedance Equation (fTl). 



3.1 A complete orthonormal system 

In ^ V. Kravchenko et al. posed a very important property of the formal 
powers, that will be the central column for our further discussions. 

Theorem 4 I5j The set of real parts of the formal powers, with coefficients 1 
and i, corresponding to (iW, valued at the boundary T of some domain fl: 



|i^eZ(")(l,0;z)|^,i^eZ(")(^,0;z)|^|" , (19) 

L J n— 

constitute a complete system for the solutions of the Dirichlet boundary value 
problem of 

This is, any boundary condition u\r can be approached asymptotically by 
virtue of the linear combination of the elements belonging to (19): 

(N N \ 

u\r - J2 ci'^ReZ(") (1, 0; z) |r - ^ c^^^ReZ^") {i, 0; z) |r = 0, 
n=0 n=0 / 

where the coefficients \cn\cn^\ are all real constants. 

L J n=0 

As a matter of fact, it is also possible to introduce an inner product for the 
set ([T]), as it was stated in [5]. In our particular case, this inner product will 
coincide with the Lebesgue integral, due to the boundary F is the perimeter 
of the unitary circle. Hence, if /(/) and g{l) are two continuous real-valued 
functions defined over F (this is I E (0, 27r]), their inner product will have the 
form ^ 

(/(0,5(0)= / ^ f{l)9m- (20) 
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Therefore, we can always obtain an orthonormal set of functions {uk}'^Q 
using the elements of (19), a property that will result most useful for our pur- 
poses. 

Specially talking about physical applications, by employing adequate numer- 
ical methods, as those described in [3] ; we will be able to approach the first 2A'' 
formal powers of the set (19), from which it will arise an orthonormal the set 
of functions {uk}'k=o^ ^ ready to approach solutions for the Dirichlet problem of 
0. 

The effectiveness of these numerical calculations has been tried in several 
works (see [5], [B] and [H]). Still, it will be very convenient for us to analyse 
here a pair of particular examples, before paying all our attention into more 
general cases. 



4 Brief study of the Dirichlet boundary value 
problem, when a is separable- variables. 

4.1 The case when a has exponential form. 

Let us consider a conductivity function of the form 

a = e^'+y. (21) 

In order to obtain at least one analytic solution for the Electrical Impedance 
Equation, we will assume that ([T]) can be separated into a decoupled pair of 
Sturni-Liouville equations: 

dy {e^'+ydyu) = 0. (22) 

It is remarkable that the modern Pseudoanalytic Function Theory has also pro- 
vided new results for this classical kind of equations. Indeed, we could even 
approach the general solution for every equation of ( |22[ ), in terms of Taylor se- 
ries in formal powers, employing the results posed in jlOj. Yet, at the moment, 
we are only searching for particular solutions. 

For this purpose, let us introduce in the first equation of (22) the auxiliary 
notation Ux = dxU. We will have that 

dxUx +Ux ^0, 

for which a particular solution is 

being f{y) an arbitrary function depending upon y. Moreover, applying the 
antiderivative with respect to x to the last expression, we can easily verify that 

u = /(y)e-^ (23) 
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is a particular solution for the first equation of (22). An identieal procedure will 
show us that 

u = g{x)e-y, (24) 
where g{x) is an arbitrary function of x, is a particular solution for the second. 



Once we combine ( 23 ) and ( 24 ) , we will immediately obtain that 

u = e-'^-y (25) 

is an exact solution of the equation ([ij . 

We have shown the full procedure because it will be useful for obtaining the 
required analytic solutions in every case we will study hereafter. 

Let us now summarise this result in a short Proposition. 



Proposition 1 Let a — e^+^. Then the function u = e ^ ^ , posed in (25), 
will he a particular solution of 

Generating a boundary condition w|r from the particular solution shown in 
the previous Proposition, is a very simple task. Remembering that F coincides 
with the perimeter of the unitary disk, we only have to consider x = cos 9 and 
y = sin 61, for G [0,27r). 

We can now employ the numerical methods detailed in [3] for obtaining a 
system of 22 formal powers, defined at the boundary F: 

|4")(l,0;z)|r,z(")(z,0;z)|r}'° , (26) 

L J n— 

that correspond to the exponential conductivity a — e^^^ . We construct each 
formal power to be composed by 1000 complex values, associated with 1000 
equidistantly distributed points on F. Then, after applying a standard Gram- 
Schmitt orthonormalising process over the real parts of ( |26[ ), we will obtain a 
system of 21 base functions {wfe}^*LQ (the apparent loss of one function can be 
explained using the expression ([8| introduced in Theorem |2j and the Defini- 
tion pi The formal power Z^'^\i,Q] z) is purely imaginary, thus its real part 
ReZW(i,0;z) =0). 

Hence, every function Uk will posses 1000 real values along F. If needed, we 
can always use the well known cubic spline interpoUing method on every Uk^ to 
posses a set of continuous functions instead of a discrete one. 

Our next step is to test the effectiveness of the set {wfclfeLo approach 
solutions of the Dirichlet problem for ([l]), imposing the analytic solution (251 



as the boundary condition u\-c- Specifically, we will evaluate ( p5| at the points 

/27r\ /27r\ f 2ti\ ( 2tx 

cose, y = sme; = 0, — , 2 — , 3 — , 20 ' 



^2iy' V21/ V21/ V 21 

and the same will be done for every function of the set {ufcj^^Q. This will 
provide a basic linear system from which we will obtain 21 constant coefficients 
{akj'l'Lo such that 

20 
fc=0 
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Table 1: Values of the coefficients corresponding to the boundary value 
problem with a of the form (21). 



ao 


ai 


02 


as 


an 


"12 


"13 


ai4 


43.708 


-13.435 


-3.713 


-1.152 


41.760 


-19.752 


2.693 


0.327 



We shaU point out tliat the number of evaluation points, located on the 
boundary F, does not necessarily have to coincide with the number k of base 
elements. We selected to do this in behalf of simplicity, but other techniques 
are available if required. The Collocation Method, successfully employed in [6], 
is an example of the last assessment. 

Remark 3 The 21 elements of the set {wfel^^g '^^'''^ obtained by applying the 
Gram-Schmitt orthonormalising process, considering first all the elements of the 
set 

|i?eZ(")(l,0;z)|r|'° , (27) 

J 71 — 

and consequently the elements of 



{i?eZ^")(i,0;z)|r} 



n=l 



(28) 



As expected, the absolute error £ will be defined according to the Lebesgue 
norm: 

1 

2 



£ = 



' .27r / 20 \ 

J (^afc?^fc-w|rj dl 



(29) 



The Figure 1(a) shows the plot of the boundary condition u|r and the plot of 



the approached solution, being the horizontal axis the perimeter of the unitary 
disk. As the reader can appreciate, it is technically impossible to detect any 
difference in plain sight between these two graphics. That is why, when this case 
occurs, we will not display the illustration, focusing our attention into some of 



the most relevant coefhcients Uk, and the absolute error £ defined in (29 1 



The semilogarithmic Figure 1(b) illustrates the absolute values of the 21 
coefficients employed for approaching the condition u|r, as well the Table 
[l] displays eight representative values of the set {uk}^^Q, keeping in mind the 
Remark |3j and rounded them to three decimal places. In other words, the first 
four values are the coefficients ak corresponding to the orthonormal functions 



upcoming from the set (26), whereas the remaining values correspond to the 



coefficients of the orthonormal functions obtained from (271 



From Figure 1(b) we could venture to say that the elements shown in Table 
[T] are among the most significant for the approach. 

The absolute error introduced in ( 29 1 , and also rounded to three decimal 
places, is £ = 2.006 x 10~®. Its integral expression was approached by the 
trapezoidal method, considering 1000 equidistant segments among the interval 
[0, 27r). This parameters will be kept for further calculations. 



10 











— Condition 
— Approach 



















































1 



2 3 4 5 6 

(a) Boundary condition and approached solution. 
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(b) Absolute values of the coefficients a^. 

Figure 1: Boundary condition and approached solution for the boundary value 
problem: Example 1. 
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Figure 2: Absolute values of the coefficients ak, corresponding to the case of 



the Lorentzian conductivity ( 30 1 



Table 2: Values of the coefficients ak corresponding to the boundary value 



problem with a of the from ( 30 ) 





ai 


a2 


as 


an 


ai2 


"13 


ai4 


40.216 


-25.102 


-28.884 


13.505 


28.106 


-13.119 


4.576 


-0.565 



4.2 The case when a has a Lorentzian form 

This example was selected because its exhibited numerical properties, that will 
be properly shown soon after. Beside, this conductivity could be situated very 
close of real Medical Imaging examples. Following identical logic steps to those 
exposed in the last section, we can pose the following statement. 



Proposition 2 Let 



Then the Junction 



a;2 + 0.1 ) + 0.1 



u = — ^+0.1(a; + 2/), 



(30) 
(31) 



will he a particular solution of 



As in the previous case, we will use (311 to construct the boundary condi- 
tion M|r, we will approach a set of orthonormal functions {1*^.}^^^, and we will 
calculate the coefhcients to approach the boundary condition. 

It is evident from the Figure 2 that the values of the significant coefficients 
ak are higher than those obtained for the exponential case, but the number of 
significant coefficients remains, somehow, the same. This can be corroborated 
watching the contents of Table [2] 

For this case, the absolute error was £ = 2.15 x 10^^. The error is consid- 
erably bigger than the one obtained in the previous subsection. Nevertheless, 
there is no need to illustrate the plots of the condition and of the reconstruction, 
since not any difference can be appreciated. In the opinion of the authors, the 
approached solution can still be considered adequate. 



12 



5 Analysis of boundary value problems when 
the conductivity is not originally a separable- 
variables function. 

Basically, all results showed above are valid if and only if the conductivity 
function cr can be expressed in terms of a separable variables function of the 



form (16). This, of course, will seldom happen in physical applications, but it 
seems to be an alternative to approach any given conductivity function in terms 
of a piecewise separable-variables function, that might be useful for practical 
applications. 

It is not clear if the proof of completeness for the set of formal powers within 
bounded domains, elegantly posed by Kravchenko et al. in [5] , can be extended 
for the case of a piecewise separable- variables function, nor it is the intention of 
this work to try to do so. Yet, the numerical results indicate that this alternative 
is worth of consideration, and the following paragraphs will intend to support 
this assertion. 

The idea for introducing a piecewise separable-variables function was first 
presented in p31 , and a basic characterization was posed in [T5] and [TH] . Here, 
the proposal will be analysed by studying six different examples. It will be 
used to approach solutions for the two cases previously studied, and it will be 
employed to analyse four cases where the conductivities are not, by definition, 
separable-variables functions. Just as performed before, exact solutions will be 
used for generating the boundary conditions to achieve in each example. 



5.1 Construction of a piecewise separable- variables con- 
ductivity function 

Consider a bounded domain f2 (in this case the unitary disk), and divide it 
into a finite number of subsections, taking care that the point to be considered 
the center of the formal powers zq (see Definition [5]) does not reside onto the 
boundary of two or more subsections. In behalf of simplicity, let us make the 
division by employing a finite set of parallel lines to the y-axis, equidistant one 
to each other, and let us locate zq = 0. 

Supposing that the values of the electrical conductivity are defined at every 
point inside the domain Q, let us trace a straight line within every subsection, 
watching out that does not intersect the bounding parallel lines of its corre- 
sponding subsection, inside fl. Indeed, such lines can simply be parallel to the 
bounding ones. 

The next step is to collect a finite set of values over every line that crosses 
the subsections. For every crossing line, the quantity of collected values must be 
large enough to warrant that an interpolling process (in our case cubic splines) 
will adequately approach all the remaining conductivity values over the line. 

Since we already assumed that every crossing line will be parallel to the 
subsection-bounding lines, and in consequence to the y-axis, all collected points 
corresponding to the same crossing line will posses the same x-coordinate. Let 
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us now propose that the conductivity inside every subsection can be represented 
according to the expression 



K 



X + K 



f{y), 



(32) 



where x denotes the x-coordinate that is common to all points along the crossing 
line, f{y) is an interpolling function that approaches the values of the conductiv- 
ity collected over the line, and K is a positive real constant such that x + K 
within the subsection. 

From this point of view, and supposing we have M subsections, the con- 
ductivity (T inside the bounded domain can be approached by means of the 
piecewise function 



My) 



cr{x,y) = 



x+Ki 
Xi+Ki 

Xm+Km ■ fMiy) ■■ X e [X^M) , X^M+1) 



■ X G [a;(i),a:(2)); 
: X e [a;(2),a;(3)); 



(33) 



Here a;(i) represents the first x-coordinate found within the domain when 
broaching the x-axis from — oo upto oo, whereas X(^m} represents the last one. 
The pairs of coordinates (a;(j), a;(j+i)), where j = 0, 1, M + 1; represent the 
common x-coordinates belonging to every pair of lines delimiting the subsec- 
tions. It is evident that the piecewise function (33) is separable- variables. 
According to the Section [Sj it immediately follows that 



Fn 



•/!(?/) 

( 



^x+K2 ' f'^^y) 



X e [x(i),X(2)); 
X e [a;(2),X(3)); 



( 



x+Km 



fMiy)) ■■ X e [X(^M),X(M+1)]- 



Whereas 



Go 



■( 



+ Ki 1 
+Ki ■ h{y) 



X+K2 . 1 
X2+K2 f2{y) 

x+K 



■ X e [a;(i),a;(2)); 
: X e [a;(2),a;(3)); 

X e [X(M),X(M+1)\- 



^Xm+Km 

For the generating pair (Fi, Gi) we will simply have 

i 



F, 
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Figure 3: Absolute values of the coefficients ak, corresponding to the case of 



the piecewise conductivity resulting from (21 1 



Table 3: Values of the coefficients corresponding to the case of the piecewise 



conductivity resulting from (21). 



ao 


ai 


a2 


"3 


an 


"12 


"13 


"14 


28.539 


-31.624 


7.445 


3.649 


44.349 


-19.430 


1.479 


0.179 



being a the piecewise separable- variables function introduced in ( 33 ) . These are 
the generating pairs (-foj Gq) and {Fi, Gi) that we will employ for the numerical 
calculations. 

5.2 The separable- variables exponential case 



Let us consider again the exponential conductivity function introduced in (21 ): 
e^+y . We will use it for approaching the piecewise-separable variables generating 
pairs {Fq,Go) and (Fi,Gi) declared in the Subsection above. 

Hereafter, we will consider M = 1001 subsections, produced by equidistantly 
parallel lines to the y-axis, and the parallel-crossing lines will be placed at the 
very middle of every subsection. Also, for every subsection, we will take 1000 
samples of conductivity, collected from an equal number of points equidistantly 
distributed over the crossing line of the subsection. Finally, we will assume 
Ki=K2 = ... = Km = 60. 

Following the methodology posed in Section [4j we can build an orthonor- 
mal system {ufe}^.^g, to approach the boundary condition u\t, raised from the 
corresponding particular solution of ([!]) m = e~'-^~y . 

The Figure 3 shows that the number of significant coefficients at has been 
increased. It seems that more than ten are among the most relevant. Neverthe- 
less, many of these values have very similar magnitudes, therefore in the Table 
[3] we will only show eight values, as is has be done before. The absolute error 
has also grown, being 8 = 3.4 x 10"'^. Still, not any difference is noticed when 
plotting together the boundary condition and the approached solution. 
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(a) Boundary condition and approached solution. 
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(b) Absolute values of the coefficients o^. 

Figure 4: Results of the calculations employing the piecewise-separable variable 
function for the Lorcntzian conductivity 



Table 4: Values of the coefficients a/, corresponding to the case of the piecewise 



conductivity resulting from (30). 
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0.057 



5.3 The separable-variables Lorentzian case 



We now consider the conductivity function (30) once more, performing the nu- 



merical calculations with all the details provided in the previous Subsection, 
and employing the piecewise separable- variables conductivity function obtained 
from (30). The boundary condition rises from the exact solution (31 1 of ([T]). 

This is the first case where it is possible to detect a clear difference between 
the plot corresponding to the boundary condition (blue line) and the approached 
solution (red line), as it is illustrated in Figure 4(a) Moreover, even the number 
of significant coefficients remains low (about four only), they do not appear 
following the same sequence that tables shown before, as we can appreciate in 
Figure |4(b)[ Please examine carefully the numbers of the subindexes contained 
in Table 4. About the absolute error, we have £ = 10.2 x 10~^. 
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Figure 5: Absolute values of the coefficients ak, corresponding to the case of 



the piecewise conductivity resulting from ( 34 1 



Table 5: Values of the coefficients ak corresponding to the case of the piecewise 



conductivity resulting from (34). 





"1 


a2 


as 


an 


ai2 


"13 


"14 


31.150 


0.001 


-1.378 


-0.001 


0.001 


17.105 
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-0.093 



5.4 The non-separable variables exponential case. 

The following four examples intend to show that the methodology for approach- 
ing piecewise separable- variables conductivity functions posed in Subsection |5.1 
can be trustful even when the conductivity function was not originally separable- 
variables. 

We will aboard this task considering first an exponential conductivity. 



Proposition 3 Let the conductivity function 



a = e 



xy 



Then, a particular solution of ^ will be 



(34) 



(35) 



Once more, we use (34) for creating a piecewise separable- variables function, 
and we impose (35) as the boundary condition. The obtained results are the 
following. 

It is not possible to notice any difference between the graphic of the boundary 
condition and the plot of the numerical approach. That is why we only show the 
semilogarithmic graphic of the coefficients ak in the Figure 5. For this case, we 
only identify two relevant coefficients, but we still show eight values in the Table 
[5] as usual. The approached absolute error for this case is £ ^ 8.948 x 10^^. 
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Figure 6: Absolute values of the coefficients ak, corresponding to the case of 



the piecewise conductivity resulting from ( 36 1 



Table 6: Values of the coefficients ak corresponding to the case of the piecewise 



conductivity resulting from (36). 
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5.5 The non-separable variables Lorentzian case. 

Proposition 4 Let the conductivity function have the form 

_ 1 

{x + yf + l 

An exact solution for the equation |7p is 

u = \-x + y. 



(36) 



(37) 



We construct the piecewise separable- variables conductivity employing ( 36 ) 



and we introduce the boundary condition using (37). For this example we 
neither find any remarkable difference between the graphics of the condition 
and the approach. The Figure 6 displays the magnitude of the coefficients ak- 
We detect only four relevant coefficients, as assented in Table [6j The absolute 
error is £ = 1.4 x lO^'^. 

5.6 The non-separable variables polynomial case. 

Proposition 5 Let us assume the conductivity function to have the form 

a = X + y + 10, (38) 

thus the function 



u = ln{x + y + 10) 



(39) 



will be a solution of 
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Figure 7: Absolute values of the coefficients ak, corresponding to the case of 



the piecewise conductivity resulting from ( 38 1 



Table 7: Values of the coefficients ak corresponding to the case of the piecewise 



conductivity resulting from (38). 
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310.378 


24.207 


-0.638 
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-22.382 
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0.068 



As usual, first it comes the approaching of the piecewise separable-variables 
conductivity based upon (38), and latter we stablish the boundary condition 
according to ( 39 ) . One more time, there is nothing to enhance from the com- 
parative graphic of the condition and the approach. The absolute values of the 
coefficients are displayed in Figure 7. Nevertheless this figure suggest that there 
are at least four coefficients playing a relevant role in the approach, we can 
venture to point out that perhaps only three are really taking a fundamental 
part. The Table 7 will show why we estimate that. We shall ask again the 
reader to observe the numbers of the subindexes. The total error for this case 
is£ = 9.8 X 10-^. 



5.7 The non-separable variables sinusoidal case. 

This case is the last one of our characterization, and we shall remark that it 
could offer interesting challenges for a more detailed examination, due to the 
possibility of numerical indetermination contained into its particular solution. 
At this moment, we will consider a basic example. 

Proposition 6 Let us consider the sinusoidal conductivity 

(7 = 1 + sinxy. (40) 

We can verify by direct substitution that the function 

«=(tan(f )+l)"\ (41) 

is a solution of (wl). 
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Figure 8: Absolute values of the coefficients ak, corresponding to the case of 



the piecewise conductivity resulting from ( 40 1 



Table 8: Values of the coefficients corresponding to the case of the piecewise 



conductivity resulting from (40). 
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-0.0881 



For the last time, and noticing the conductivity (40) does not vanish inside 
the unitary disk, we will employ it for approaching the piecewise conductivity 
function, and (41 1 for constructing the boundary condition. The Figure 8 il- 
lustrates the absolute values of a^, whereas the Table 8 shows us that, for this 
example, only two coefficients take an important part into the approach. The 
total error is f = 7.694 x 10"''. 



6 Conclusions 

Based upon the results showed in the previous sections, we consider that the 
basic characterization of the new numerical method for approaching solutions 
of the Electrical Impedance Equation is acceptable, since it has provided the 
proper information for performing new experiments on a wider class of con- 
ductivity functions that, by definition, are not separable- variables, as well as 
for starting the analysis of geometrical conductivity distributions, whose exact 
mathematical representation are, in general, unknown. 

This second class of conductivity distributions is among the most interesting 
for several branches of Experimental Physics, as it is the Medical Imaging. 
But it is clear that the characterization of the new method over this class of 
conductivities, can only be adequately achieved by comparing its numerical 
results with those provided by some of the finest variations of the Finite Element 
Method, that have proved to be the best tool to approach solutions for Dirichlet 
boundary value problems in the plain. 

The selection of the proper examples for these tests could well be worth of 
a complete paper. Anyway, it is absolutely necessary to perform those compar- 
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atives. 

We must also emphasize that we have studied only one class of piecewise 
separable-variables conductivity approach. Indeed, to try again the posed ex- 
amples, and as much of their variations as it is possible, using other ideas for 
obtaining separable-variables expressions at every subsection, will enrich our 
understanding of this new technique from the mathematical point of view. 

On the same direction, it is also necessary to formally extend the proof 
of completeness provided in [5], if such was possible, for piecewise separable- 
variables functions. About this, the authors do not detect at the present any 
immediate way to start accomplishing the task. 

And we still have to take into account all material related to the Compu- 
tational Complexity of the new posed numerical methods, that was completely 
out of the scope of the current paper. Yet, this path could be somehow clearer 
to start exploring, if we keep in mind that our ultimate objective is to apply the 
methods for Medical Imaging problems, so it will be obligatory to include the 
computational cost and time as two parameters for evaluating the technical effi- 
ciency of the new methods. This because the known algorithms for approaching 
solutions for the inverse problem, are based on the recursive resolution of the 
direct Dirichlet boundary value problem. 

Still, we are certain to be searching in the correct direction, because the tech- 
niques provided by the Pseudoanalytic Function Theory, classical and 
modern [3], are proving to be specially powerful in Mathematical Physics, e.g., 
these techniques could well prove useful in an assortment relevant applications, 
as it is the Fourier transform of the two-dimensional Fokker-Planck equation, 
studied in [8!, as well as some certain special cases of the diffusion equation 
employed into hydrodynamic fluctuations |2- 
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